A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

08 (01.01.2; 02.02.1; 13.07.2; 08.02.3) 



ASTRONOMY 

AND 

ASTROPHYSICS 

1.2.2008 



OS 



O 

CO 



CN' 

>■ 
Os\ 
CO . 
O- 

o- 

OV 
5' 

Or 
i ' 

o 



S 



The 1 - 12 Hz QPOs and dips in GRS 1915+105: tracers of 
Keplerian and viscous time scales? 

S. Trudolyubov 1 , E. Churazov 2 ' 1 , M. Gilfanov 21 

1 Space Research Institute, Russian Academy of Sciences Profsoyuznaya 84/32, 117810 Moscow, Russia 
Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85740 Garching bei Munchen, Germany 



Abstract. We analyzed 9 RXTE/PCA observations of 
GRS 1915+105 in the flaring state, when the hardness of 
the source spectrum was changing on the time scales from 
few seconds up to 1000 seconds. The quasi periodic oscilla- 
tions (QPOs) with the frequency varying between ~ 1 and 
12 Hz are associated with the episodes of harder source 
spectrum. In each observation we found tight correlation 
between the duration of the hard episode and the charac- 
teristic QPO frequency. For a half of the observations this 
correlation matches the relation between the viscous time 
scale and the Keplerian frequency, when both quantities 
are evaluated for various radii in the radiation pressure 
dominated accretion disk. Assuming that the QPO fre- 
quency is proportional to the Keplerian frequency at the 
boundary between an optically thick accretion disk and a 
hot comptonization region, the changes of the QPO fre- 
quency can then be understood as due to variations of this 
boundary position on the viscous time scales. 



Table 1. The list of RXTE/PCA observations of GRS 
1915+105 used for the analysis. 



Obs.ID 


Date, UT 


Start, UT 


Exp. a , s 


10408-01-01-01 


06/04/96 


05:39 


5068 


10408-01-10-00 


26/05/96 


17:12 


6128 


10408-01-38-00 


07/10/96 


05:44 


9766 


10408-01-44-00 


25/10/96 


11:52 


8250 


20402-01-01-00 


07/11/96 


05:42 


6948 


20402-01-28-00 


18/05/97 


16:19 


7349 


20402-01-33-00 


18/06/97 


12:58 


6093 


20402-01-35-00 


07/07/97 


14:53 


5718 


20402-01-59-00 


17/12/97 


02:10 


8795 



Deadtime corrected value of the PCA exposure 



2. Observations and data analysis 
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1. Introduction 

The X-ray source GRS 1915+105, one of the sev- 
eral galactic objects producing observable superlumi- 
nal jets (Mirabel & Rodriguez 1994), was discovered by 



GRAN AT observatory as a transient in 1992 (Castro- 
Tirado, Brandt & Lund 1992). Long-term monitoring of 
GRS 1915+105 with RXTE revealed repetitive character 
of the source variability patterns. Basing on the results 
of spectral and timing analysis, bursting behavior of the 
source can be roughly reduced to the sequence of varying 
'hard' and 'soft' states qualitatively distinguished by their 
spectral and temporal properties (Belloni et al. 1997a, 6; 
Markwardt et al. 1999; Muno et al. 1999). 



The observations used for the analysis are listed in Table 
1. For the timing analysis in the 2 — 30 keV energy range 
stars: individual: thc RXTE/PCA (Bradt, Swank & Rothschild 1993) data 
in the 'binned' and 'event' modes containing X-ray events 
below and above 13 keV respectively were used. We gener- 
ated power density spectra (PDS) in thc 0.25 — 150 Hz frc- 
qucncy range accumulated in the 4 s time intervals (Leahy 
et al. 1983). The PDS were linearly rebinned into 0.25 Hz 
bins. The white- noise level due to Poissonian statistics cor- 
rected for the dead-time effects was subtracted. For the 
determination of the QPO centroid frequencies the PDS 
were fitted to the analytic model consisting of a Lorentzian 
profile and a power law continuum in the 0.5 — 15 Hz fre- 
quency range (typical error of the derived value of the 
QPO frequency is - 0.2 - 0.3 Hz). 
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3. Results 

Two typical examples of GRS 1915+105 X-ray flux his- 
tories in the 2 — 30 keV band are presented in Fig. [j] 
along with hardness ratios (13 — 30 keV to 2 — 13 keV) 
and dynamic PDS integrated over 4 s time intervals. The 
source behavior in the flaring state is characterized by fre- 
quently occurring episodes of hard energy spectrum last- 
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GRS 1915+105; 07/10/1996; START TIME = 06:36 
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Fig. 1. Part of the light curves (upper panels), corresponding hardness ratios (13 — 30 keV)/(2 — 13 keV) (middle panels) and 
dynamic PDS (lower panels) of GRS 1915+105 for the Oct. 7, 1996 (left panels) and Oct. 25, 1996 (right panels) observations 
(2 — 30 keV energy band, PCA data). The QPO peak appears as 'U'-shaped dark band in the lower panels. The light curves were 
not corrected for the instrument dead time which was 5 — 20%. The overall count rate corresponds to 5 Proportional Counter 
Units of the PCA detector. 
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Fig. 2. The relation between the duration of a hard episode and a corresponding minimal QPO frequency for the flaring state 
observations listed in Table 1 (first group - left panel; second group - right panel). The dependences £ v i sc oc f K of the viscous 
time scale upon the Keplerian frequency at the inner edge of the radiation pressure dominated disk with the mass accretion rate 
m ~ 0.11(a/0.1)" 1 ' 2 (m/33)~ 2 ' 3 and m ~ 0.17(a/0.1) -1 ' 2 (m/33) -2 ' 3 are shown in the left and right panels by long-dashed and 
dotted lines respectively (see text). 
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ing typically 10 — 1000 s. The prominent, relatively narrow 
variable-frequency 1 — 12 Hz QPO peak (dark 'U'-shaped 
band in the lower panels in Fig. 1) is a generic feature of 
the source power density spectrum during hard episodes 
(Markwardt et al. 1999; Muno et al. 1999). 

Using the data from our sample of flaring state ob- 
servations (Table 1), we studied the relation between the 
duration of a hard episode, th ari i (determined using the 
hardness ratio) and the lowest value of the QPO frequency 
reached (determined from the fitting of the PDS to the an- 
alytical model). In Fig. the duration of a hard episode 
is plotted as a function of the minimal QPO frequency. 
For all observations used in the analysis there is a strong 
anticorrelation between these two quantities. The depen- 
dence has a nearly power law form for durations of a hard 
episode shorter than ~ 100 — 200 s and probably flattens 
for longer events. According to the value of the slope F of 
this power law dependence the observations can be sub- 
divided into two groups: T ~ 2.1 - 2.4 for the 06/04/96, 
07/10/96, 18/06/97, 07/07/97 and 17/12/97 observations 
(the first group, left panel in Fig. g) and T > 3.0 for the 
26/05/96, 25/10/96, 07/11/96 and 18/05/97 observations 
(the second group, right panel in Fig. |2|). The analysis 
of the flux histories showed that these two groups have 
qualitatively different lightcurves: for the first group the 
transition from 'hard' to 'soft' state is characterized by the 
growth of the 2 — 30 keV flux (Figure [fl upper left panel) , 
while for the second group this transition corresponds to 
the decrease of the X-ray flux (Figure EL upper right panel). 
In the subsequent paper (Trudolyubov et al. 19996) it will 
be shown that these two groups are also distinguished by 
the properties of their energy spectra: for the first group 
the soft thermal component makes a significant contribu- 
tion to the overall luminosity in the 3 — 20 keV energy 
range, for the second group the energy spectrum is domi- 
nated by the hard component (the contribution of the soft 
thermal component to the total luminosity is < 25%). 

4. Discussion 



lar to that expected if the duration of the hard episode 
is proportional to the viscous time scale of the radiation 
pressure dominated disk (Belloni et al. 1997&). Given (i) 
and (ii) some sort of correlation between the QPO fre- 
quency and the duration of the hard episode is naturally 
expected. 

To explain the observational properties of black hole 
and neutron star binaries in the hard state a num- 
ber of models involving the hot comptonization re- 
gion near the compact object surrounded by the opti- 
cally thick accretion disk was proposed. It is often as- 
sumed that the QPO phenomenon is caused by inter- 
action between these two distinct parts of the accre- 
tion flow occurring on the local dynamical time scale 
at the boundary of these regions ( Molteni et al. 1996; 
Titarchuk, Lapidus fc Muslimov 1998 ). In the following 



analysis we will also assume that the QPO frequency 
is proportional to the Keplerian frequency at the inner 
boundary of the accretion disk. 
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The study of the spectral and temporal evolution of GRS 
1915+105 during the periods of flaring activity has shown 
that bursting behavior of the source can be reduced to 
a sequence of varying hard and soft episodes. The hard 
episodes are distinguished by the presence of the promi- 
nent QPO peak in the source power density spectrum. 
Combined time-resolved spectral and timing analyses re- 
vealed (i) a strong correlation between the QPO frequency 
and parameters of the soft component of the energy spec- 
trum (presumably emitted by the optically thick part of 
the disk) on a wide range of time scales (Trudolyubov 
et al. 1999a; Markwardt et al. 1999; Muno et al. 1999); 
(ii) the dependence of the duration of the hard episode 
upon the corresponding characteristic radius of the opti- 
cally thick disk (derived through the spectral fitting with 
multicolor disk black body model) thard oc r 7 ' 2 is simi- 



Fig. 3. Dynamic power density spectra for two represen- 
tative observations (18/06/1997 - first group (upper panel, 
25/10/1996 - second group lower panel). The expected time 
dependence of the QPO frequency implied by viscous evolu- 
tion of the radiation pressure dominated part of the accretion 
disk is presented by a white solid lines superimposed on the 
observed QPO track. 



In the standard accretion disk theory (Shakura & Sun- 
yaev 1973) the viscous time of the radiation pressure dom- 
inated disk is ivisc ~ 1-2 x 10 _5 a _1 mrh _2 r 7 / 2 s, where 
a - viscosity parameter, m - mass of the compact object 
in solar masses, m - disk accretion rate in units of criti- 
cal Eddington rate, r - distance to the compact object in 
units of 3 gravitational radii. Introducing the local Keple- 
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rian frequency, /k 



2200 m,- 1 ^ 3 / 2 Hz, we obtain: 



U 



740a 



_1 m 4,/3 to 



-2,-7/3 



(1) 



-7/3 



This dependence, t v i sc oc / K , reproduces well the 
relation between the duration of a hard episode and the 
associated minimal QPO frequency for the first group of 
observations (Fig. @, left panel). Moreover, for several ob- 
servations from this group the observed relation between 
the duration of a hard episode and the maximal value of 
the characteristic inner radius of the disk, determined via 
spectral fitting, matches remarkably well the expected de- 
pendence of the viscous time scale on the radius for the 
radiation pressure dominated disk (Bclloni et al. 19976). 

Assuming that the QPO frequency is proportional to 
the Keplerian frequency (/qpo = /k/0 and duration 
of the hard episode is proportional to the viscous time 
(^hard = tv\sc/k), one can estimate the required mass ac- 
cretion rate m using equation (1): 



a \ - 1 / 2 ( m \ ~ 2 / 3 ( k ih. 



01/ V33/ 



hard 



100s 



-1/2 



(z/ Q p r 7/6 (2) 



In Eq. (2) compact object mass was normalized to the 
value of 33 M , implied by the interpretation of the 
67 Hz QPO feature as a signature of a Keplerian os- 
cillations near the last marginally stable orbit around a 
Schwarzschild black hole (Morgan et al. 1997). The in- 
ferred values of the accretion rate for the observations 
from the first group are in the range to ~ (0.11 — 
0.17)(a/0.1)- 1 / 2 (m/33)- 2 / 3 fc- 1 / 2 r 7 / 6 (Fig. §. This con- 
clusion is supported by the results of the spectral analysis 
(Trudolyubov et al. 1999&) where the bolometric luminos- 
ity of the soft spectral component corresponding to a given 
value of QPO frequency was used as a measure of the mass 
accretion rate. 

If the QPO frequency is associated with the selected 
region in the accretion disk (e.g. inner edge of the optically 
thick disk) which is moving radially on the viscous time 
scale, then we expect correlation of the QPO frequency 
and the rate of the frequency change. Assuming that rise of 
the QPO frequency is caused by the inward motion of the 
inner edge of the radiation pressure dominated accretion 
disk one can write: 



dt 



dr 



v T (r) 



■ 2 7/3 j-10/3 

am to ' / QP ' , 



(3) 



where v r (r) - radial velocity of the matter in the disk at 
radius r (Shakura & Sunyaev 1973). Integrating Eq. (3), 
we obtain: 



/qpo (*o) ~ /qpo (*) - Mt ~ *o), 



'QPO 



(4) 



where A ~ am 2 m 7 ' 3 . In Fig. |3| this dependence is shown 
in comparison with the observed QPO evolution for two 
observations from the first (18/06/1997) and the second 



(25/10/1996) groups. For the first group both initial de- 
crease and following rise of the QPO frequency during the 
hard episodes are generally described by Eq. (4) (Fig. [| 
upper panel) but with different values of coefficient A. For 
the rise phase the value of coefficient A remains practi- 
cally the same on a time scale of an individual observa- 
tion (~ 10 4 s). For the second group of observations only 
rise of the QPO frequency is generally described by Eq. 
(4), while the initial decay often has a more complicated 
structure (Fig. ||, lower panel). In addition, an extended 
plateau in the time dependence of the QPO frequency near 
its minimum is sometimes present. This might explain the 
deviation of the dependence of thard on the QPO frequency 
from the thard oc /q P o l aw f° r the second group of obser- 
vations. 



5. Conclusions 

We analyzed the temporal behavior of GRS 1915+105 us- 
ing the set of 9 RXTE/PCA observations representing typ- 
ical patterns of the source X-ray evolution during the pe- 
riods of flaring activity. 

1. For all individual observations the tight correlation be- 
tween the duration of a hard episode and the minimal 
QPO frequency was found. 

2. For ~half of the observations this correlation is satis- 
factory described by the standard radiation pressure dom- 
inated disk theory assuming that the QPO frequency and 
the duration of a hard episode are proportional to the 
Keplerian frequency and the viscous time on the inner 
edge of the disk respectively. The derived values of the 
mass accretion rate for these observations are in the range 
to ~ (0.1 - 0.2)(a/0.1)- 1 / 2 (TO/33)- 2 / 3 . Not all observa- 
tions fit this dependence. In fact for the large fraction the 

— 7/3 

relation thard oc /qpo breaks. It should be noted, that 
for these observations the correlation thard oc r 1 1 2 (Bclloni 
et al. 1997&) also fails. 

3. The temporal evolution of the QPO frequency during 
the rise phase of the hard episodes can be explained as 
a change of the Keplerian frequency at the inner edge of 
optically thick accretion disk caused by its viscous motion. 
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